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Abstract 

In this paper we prove the existence and uniqueness of the solution of a non-stationary 
problem that modelizes the behaviour of the concentrations and the temperature of gases 
going through a cylindrical passage of an automotive catalytic converter. This problem 
couples parabolic partial differential equations in a domain with one parabolic partial 
differential equation and some ordinary differential equations on a part of its boundary. 


1 Introduction 


The starting point of this work is a non-stationary model of catalytic converter with cylindrical 
passages due to Ryan, Becke and Zygourakis [12] that is an extension to the one established by 
Oh and Cavendish [10,11]. It describes the spatial and temporal evolutions of the concentrations 
of N — 1 different chemical species in gazeous phase going through a cylinder and that of the 
temperature in the cylinder and on its boundary. A stationary problem close to our has been 
studied by Carasso [2]. Because of its internal symmetry the cylinder may be reduced to the 
plane domain O = ]0,1[ x ]0,1[ the boundary of which is E = {1} x ]0,1[. For i G {1, ...,N — 1}, 
the concentrations (resp. the temperature) inside the cylinder O are named Cif (resp. C^f) 
and the concentrations (resp. the temperature) on the boundary E are named Cis (resp. C^s)- 
The problem is written in a normalized way as 


dC, 


if 


(r, z, t) = 


I3if 


d ( dC, 


dz 

ac, a^Cn 

( 2 , tj - ^ [Z, t) 


if 


dt 


dz'^ 


r (1 — r^) dr 
^^if n 


dr 


ir,z,t ), 


( 1 ) 


(C*ls) • • • ) C*Afs) (^) ^) 
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for i G N}, 6i G {—1,1} and = 1 if i = N, and = 0 if not, moreover the constants 

A/) 7is cind ^TVs are nonnegative. The initial and boundary conditions are 

^ dC 

Cif{r,0,t) = Cifo{r), -^{0,z,t) = 0 , 

< Cif{l,z,t) = Cis{z,t), Cis{z,0) = Ciso{z), ( 2 ) 

eNs^il,t) = 0 , = 0 . 

The functions r^, i G N}, are supposed to verify the following assumptions 

(HI) They are Lipschitz continuous with the same constant k. 

(H2) For all (xi,..., xn) in we have r* (xi,..., xn) > 0. 

(H3) If one of the Xj, j G {1,..., A^}, is equal to 0 then r* (^^i, • • •, 0) • • •) ^n) = 0. 


Remark 1 The physical meaning of these assumptions are 

1. The function r* is the speed of creation or consumption of the i — th chemical species 
through all the reactions in which the species appears. 

2. (H2) means that the i-th chemical species is either created (6i = 1) or consumed (6i = — f). 

3. (H3) say that if there is no more of one chemical species necessary for the reaction with 
the i-th species, then the corresponding reaction stops. 


Integrating the Tth equation of Oi with respect to r after multiplying by r (1 — r^) allows us 
to rewrite the system o as 


r ^ 

dz 

< dCis 


A/ 

d 


r (1 — 

) dr 

\ dr 


dC.f 

r{l-r 

A/io 

dz 

-f-diTi (Cis,... 

, Cns) ■ 


( 3 ) 


It is this problem (ED with the initials and boundary conditions 0 that we are studying from 
now. 


2 Existence of the solution 

To prove the existence of a solution of ED the main idea is to perfom what we called a “parabolic 
regularization” of (02 in order to avoid working with both pdes and odes and hrst show the 
existence of a solution of the new regularized problem. Then we deduce the existence of a 
solution of the initial problem using a limiting process. 
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2.1 Regularization of the problem 

We regularized Q by adding with 9is > 0 and i G {1,..., iV — 1}, in the equations 


®2 on the boundary to obtain 

dC.f 
d 

da 


-e,_ 


dz 

d^a 




r (1 — r^) dr \ dr J ^ 

r (1 - r^) dr + diVi , Cns) 


lis f^dCif 2 


dt dz'^ A/Jo dz 

with the initials and boundary conditions m and 
BC BC- 

(1, t) = 0 = (0, t), w e {1,..., iv -1}. 

Then we prove the existence of a solution of the problem 0, using the mappings 

T $ 

Cis Cif and Qf Qs 

Lemma 2 There exists a nonnegative constant c (T) such that 


N N „T rl 

T {c.,fe,T)dz + y2^t, / 

i=ao Jo Jo 


dq 

dz 


dzdt < c (T). 


Proof. From (Q, we deduce after summation on i and using (HI) and (H3) that 




dq 

dz 


dzdt 


< a{T) + sup ( ^ ) kNj2 


lis 


N rT rl 


i=l 



N „i 


_i JO JO 
N 


(Cjs)^ dzdt. 

A/ r 


withalTj^Wf {Q«y-{r)r{l-r-)dr + q^f (O^o) 

^ i=l -^0 ^ Jis Jo 


{z) dz, which leads to 


N A 


N rT pi 


A/ (C.,f(z,T)dz<h(T) + dVl [ (c, 

Jo .•_1 Jo Jo 


dzdt 


i=lJ0 JO 


( 4 ) 


( 5 ) 


with 6(T) = 2a(T)/infi (^) and d = 2supi (^) kN/infi (^) . Using Gronwalhs lemma 

y / y J y i 

leads to the result. ■ 


2.2 Existence in the cylinder 

Assume that Cis, N} are known on the boundary (mapping T) and let 

Uf = * {Cif,..., C^f) , Ufo = * {Cifo, .. ., Cisifo) 5 

Us = * (C*is, ..., Cns) , (Jf = diag (A/, • • • , (^Nf) ■, 
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to obtain the following problem in the cylinder 


duf 

dz 


-Pf 


1 d f duf 


= 0 , 


r (1 — r^) dr \ dr 

Uf{r,0,t) = Ufo{r), 

Uf{l,Z,t) = Us{z,t), 

^{0,z,t) = 0. 

Performing the change of fonction Wf {r,z,t) = Uf (r, 2 ;, t) — Ug {z, t) we have the problem 


( 6 ) 


dm, 

dz 


Pf 


1 d f dwf 


dug 


r (1 — r^) dr \ dr J dz’’ 

w/(r,0,t) = Wo(r,t), 

Wf{l,z,t) = 0 , 

dwf 


(7) 


dr 


{0,z,t) = 0, 


with Wo (r,t) =Ufo (r) - Ug (0,t) =Ufo (r) -ujo (1). 

Definition 3 1. We set 


7^^(i_r2) (0,1) = I r M (r) A/r^TT^rT) e (0,1)| , 


(0,1) = {m I r M (r) a/t e (0,1)} . 

-^r(i-r 2 ) 7) Ll{0,l)) is a Hilbert space for the norm associated to the scalar 

product 

(M,n )2 ,,(i_^ 2 ) = [ uvr [l - r'^) dr (resp. {u,v)^^^= f uvrdr). 


2. Let 


V(l-r2) (0) 7)) 


^ , du 


dr 


e {Ll ( 0 , 1 )) 


N 


Wr = ^ue [Li 

Wro = {ueWr |m( 1) = 0}, iP,(T) = {ue L^]0,l[x]0,T[-,Wr)} 
and be the dual space ofWro- 

3. The spaces Wr, Wro and Wr (T) are equipped with the norms 


u 


Wr 


U 



h 

du 

L 

dr 


du 


dr 


rdr, 


rdrdzdt, 


where ||.|| is the Euclidian norm ofM.^. 
We establish some properties of these spaces 
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Lemma 4 1. We have the continuous embeddings 


(0,1) C Ll (0,1) C (0,1). 

2. For all u in Wro, we have 



(8) 


/ N Af 

Then the injection of Wro in (0)1)] continuous and compact and the em¬ 

bedding is dense. 


3. Wro is a Hilbert space for the norm above. 

Proof. 1. As Vr G [0,1], 0 < r (1 — r^) < r < 1, we deduce the result from Kufner [9] . 
2. For all a such that 0 < a < 1 we have as m( 1) = 0 that 


u(a)|| 




du 

dr 



Using the Cauchy-Schwarz inequalitie, Integrating and multiplying by a(l — a^) leads to the 

/ N w 

result. We show the compactness of the injection of Wro in ( (0,1) j and the density 

of the embedding as in Dautray-Lions [4]. For more details, see Hoernel [6]. 

3. Immediate consequence of 2. ■ 


Definition 5 Let Ufo G (0,1) j and Ug G ^0, T; (0, l))^j . A function Wf 

is called a weak solution of ((7j) ^f and only if e lU, (T), ^ G (]0,1[ x ]0, T[; 

Wf (r, 0, t) = wq (r, t), and if for every test-function (p G (]0,1[ x ]0, T[; Wro) we have 



T 


■ ) r (1 — r^) drdzdt + 

(fr (1 — r^) dr ) dzdt 



0 JO JO 


dwf dip , , , 

■ - 7 :— r dr dzdt 
dr dr ' 


(9) 


To prove the existence of a weak solution of dZD we dehne the second order elliptic operator 
L : Wro ^ W;, by 

f dw dv\ 

Ww,vE IF,o : j rdr 


N 


and the operator T from (0,1) j in itself giving for each function g taken in (^F^|.^_^ 2 ) (0,1) 

the unique element T [g) in Wro that is the weak solution of F (T [g)) = gr {1 — r^) by 




dTjg) 

dr 


dip 

dr 


^ rdr = J {g ■ ip)r (1 — r^) dr, 'iip G IWo- 


( 10 ) 
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Lemma 6 1. The operator T is 

- Self-adjoint: V/, g G (0,1) j : (T (/), = (/, T 

- Positive: V/ G (0,1) j : (T (/), /)^(i_^ 2 ) > 0 

- Non-degenerated: V/ G (0,1) j : (T (/), /)^(i_^ 2 ) = 0 ^ / = 0. 


Moreover there exists a Hilbert basis {uj^}j of both 
eigenfunctions of T and L. 


\Ll{i-r 2 ) (0,1)j and Wro made of 


Proof. 1. Immediate. 

/ X AT 

2. As the injection of HAo in (-h^(i_r 2 ) (0,1) j is compact, T is a compact operator from 


-^r(l-r2) (0) 1) 


N 


in itself. Taking into acount the properties of T, the Theorem VI. 11 from 

/ \ N 


Brezis [1, page 97] implies the existence of a Hilbert basis of (0, l)j made of 

eigenfunctions of T and of L. The Remark 29 from Brezis [1, page 193] implies that this basis 
is also a Hilbert basis of H7o- ■ 

We show the following existence result 


Proposition 7 Let uq and Us as in Definition^ Then, there exists at least one weak solution 
Wf o/dH). 

Proof. Let be a basis of HVo made of eigenfunctions of L with eigenvalues \j 



= A, 


{uj^ • (p) r (l — r^) dr, 


V(p G HAo- 


For all (p in HVo, there exists a sequence such that (p = ■ AVe write Pm {g>) = 

Galerkin approximation of order m of (p. Properties of implies that the 

operator Pm is a continuous orthogonal projector from H7o (resp. 

{Ll(^i-r 2 ) (0,1) j ) to the subspace of Wro (resp. of ( T^(i_^ 2 ) (0,1) j ) generated by a;\ ... 
Moreover 


II ( 7^)11 VKro — llT’llwrO ’ ll-^"i (T’) Il2,r(l-r2) — 117^112, r(l-r2) 5 V 79 G ITrO. 

We establish the following estimates 

,t)r{l — r‘^)drdt < C{T), 

2 

{r, z,t) rdrdzdt < C{T), 

< C{T), 



(11) 


L2(]0,l[x]0,T[;W;o) 
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for almost every s in ]0,1[, where C{T) is a nonnegative constant that only depends of the 
problem m datas and T, but independant of m. Then, there exists a subsequence of {Pm {wf))^ 
still noted {Pm {wf))^ such that 


Wf, weakly in (^]0,l[x ]0,T[; 1)) ^ , 

weakly in (^]0,1[ x ]0, T[; (L2(0, 1))^) , 
weakly in (]0,1[ x ]0, T[; W'q) . 

Consequently {Pm {wf))m weakly converges to Wf in (]0,1[ x ]0, T[; Wro) and Wf verihes the 
variationnal formulation m- ■ 

Remark 8 Ifufo belongs to Wr then E Li^ ^]0,1[ x ]0,T[; (0, l)j ^ 

Remark 9 From (]0,1[ x ]0,T[; IT^o) C ^]0,1[ x ]0,T[; (0, l)j ^ , 

^]0,1[ X ]0, T[; (0,1) j ^ C (]0,1[ x ]0, T[; MTq) and the fact that we have 

Wf E L^(]0, l[x ]0, T[; hTj-o) and E L^(]0, l[x ]0, T[; using Proposition 23.23 from 
Zeidler [13, page 422] we obtain that wj E C ^]0,1[ x ]0,T[; (0,1) j ^ which allows 

us to interpret the initial condition. 


Pm {Wf) -f 

m—>+oo 

dPm {Wf) 

dr m—>+oo 

dPm {Wf) 


dz 


m^+oo 


Corollary 10 The problem admits a weak solution Uf E Wr (T) such that 

(§7 ■ •") ’■ *** ^ I II = 

for allT > 0 and all ip E (]0,1[ x ]0, T[; Wro) and Uf verify Uf {r,0,t) = Ufo (r). 



0 , 


Proof. Let us notice that Uf = Wf + Ug is a. weak solution of (0) belonging to the right space 
and verifying the initial condition at z = 0. ■ 

Remark 11 Since ^ G (]0, 1[ X ]0,T[; Wro) and the function r i—^ r (1 — r^) is equal to 
zero at r = 1, we can affirm that for all function g : {z,t) g {z,t) in {Lf (]0,1[ x ]0,T[))^, 
the function (r, z,t) g {z, t) r {1 — r^) is in the space (]0,1[ x ]0, T[; IT^o) the following 
duality bracket is well defined 

^,9r {l-r^)\ 

/ {L2(]0,l[x]0,T[;m;o),L2(]o,l[x]0,T[;m,o)) 
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2.3 Existence on the boundary 

QC' 

We use the mapping <h. Suppose given the for z G {1,..., N} and let 


UsO 

r 

r 


—* (ClsO) • • • ) Cnso) ) 

=* (ri,..., rjv), 


= diag 



INs 
’ ^Nf 


6 = diag (^i, • • • ,6n) , 

9s = diag {9is,-■ ■ ,9 ns) , 

=: diag{Ti,... ,Tn) , 


in order to obtain the system 




"^duf 


Let 


dt dz^ “‘v-A/ - 

Us {Z, 0) = UsQ (z), 


W, = \ue (ryo.l))" I (lUo,i))"|., W,{T) = {u e U{(,,T-,W,)} 


r (1 — r^) dr, 


( 12 ) 


dz 


and be the dual space of W^. These spaces are equipped with the following norms 


\u\ 




'0 


||m||^ dz + 

T /•! 



du 

/„ 

dz 


dz, 


\u\ 


W.{T) 



T rl 


\u\\ dzdt + 


0 JO 



0 JO 


du 


dz 


dzdt. 


Definition 12 Assume ^ G (]0,T[ x ]0,1[; W^q) and Uso G (0,1))^ . A function Us is 
called a weak solution of m if and only if we have Us G Wz (T), Us {z, 0) = Uso (z), and if, 
for all xjj in (0, T; Wz), we have 


ff 


■ f) ) dzdt + 



' dus 
{6r{us) ■ 'ip) dzdt 



0 Jo 


' " 'efff . f) i-Mt 

'0 JO \ J 

'T pl / Q \ 

(1 — r^) \ dzdt, 

'0 JO \ / (W'^o,Wro) 



(13) 


for T <T, cf. Rem.arkrrn. 


Proposition 13 Letuf anduso be o,s in the Dehnition^VB. Then, there exists at least one weak 
solution Us of (j12jl . 

Proof. We use some auxilliary linearized equation and some fixed point as in Chipot [3] or 
Henry [5]. For more details, see Hoernel [6]. ■ 
























2.4 Existence 

We begin by showing some properties of and ib 
Proposition 14 The mapping 

^ / WriT) W^T) 

Uf Us 

which for every Uf E Wr (T) gives the weak solution Ug G Wz (T) of (IT^ verifies 

b (“/-“/) 




W,{T) 



0 JO 


dz 


dzdt 


(14) 




where a and b are two nonnegative constants independant on T. 

Proof. Let Ug = u\ — be the difference of two solutions of m with same initial conditions 
at t = 0 associated to «]• and uj- in Wr (T). Ug is a weak solution of 


— Og ^ ^ = — / T-^r {1 — r^) dr + 6 {r{ul) — r (uD) , 


dt dz'^ j 0 

= 0 , 


(15) 


with Uf = — u'j. Multiplying (USDi by Ug, integrating on ]0,1[ x ]0, r[ with 0 < r < T and 

using (HI) leads to 

2 


\[ \\Ugf{z,T)dz + 


'0 


ff 


mdUg 


f*T rl pT pi 

< k I ||t/s||^ dzdt + c 




dz 

.dUf 


dzdt 


dz 


II 17^11 dzdt 


Ko 


>0 Jo Jo Jo 

with c > 0. Using Young’s inequality with e > 0 and Gronwall’s lemma we have 


T rl 



0 JO 


dUg 


dz 


dzdt < 


cee 


{c+ 2 k)T/e pT p 1 


2 inf, 9. 



Jo Jo 


dU 


and taking e: = (c + 2fc)T we hnally obtain 

ll^sllw^(T) — + 2^) (e — 1) + 


c(c + 2k)e 
2 inf, 6^ 


T 


dz 


T pi 


dzdt 


w' 


/Jo Jo 



.dUf 


dz 


dzdt. 


WL, 


Proposition 15 The mapping 


^ : 


1U,(T) ^ Wr{T) 

Ug I—*• Uf 


which for every Ug G Wz (T) gives the weak solution Uf of ^ is such that 

11'^ i'^s) ~ ^ (■^s) |lvi/r(T) — ll^s ~ ^«llwz(T) 
for a nonnegative constant c that only depends on f3if. 


(16) 
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Proof. Let t// = wj- — be the difference of two solutions of corresponding to u], and 
ul of Wz (T), with same initial conditions at z = 0 and set Ug = u] — uj. and Wf {r,z,t) = 
Uf (r, z, t) — Ug {z, t) to obtain 


r dWf 

dz 




r (1 — r^) dr 


d f^dWf 
dr 


dUg 
' dz ’ 


iy/(r,0,t) = 0, 
Wf{l,z,t) = 0, 
dWf , 

-g^(0,j,i) = 0. 


(17) 


Multiplying (ED by r(l — r‘^)Wf, integrating on ]0,1[ x ]0, s[ x ]0, r[ with s G ]0,1[ and r G 
]0, T[ we obtain using Young’s inequality for e > 0 and r (1 — r^) dr = 1/4 that 


r /»! 



T /‘S /»1 


||1L7||^ (r, s, t) r(l — r‘^)drdt + 2 


0 JO 



e 

< - 
- 4 


T nl 



0 ^0 


dUg 


dz 


T ns nl 


0 ^0 ^0 




i/2dWj 


f Qj. 


rdrdzdt 


dzdt H— 



||1L/||^ r(l — r‘^)rdrdzdt 


0 JO JO 


and we deduce the following inequalities from Gronwall’s lemma 


r ns nl 



'0 ^0 ^0 

ns nl 


||1L/||^ r (1 — r^) — l) 


dz 




0 JO JO 


P 


i/2dWf 


f Qj. 


rdrdzdt < 


r nl 


0 ^0 


dzdt, 



0 JO 


dUg 


dz 


dzdt. 


Taking e = s we notice that because Ug doesn’t depend on r to have 


dr 


r t f 

d{^{ul)-^{u^g)) 

" e r 

d - ul) 

Jo Jo Jo 

dr 

1 Lit LL/JLLL ^ n / ^ \ 1 1 

- 8inf,(A/)yo Jo 

dz 


dzdt. 


From the fact that Uf = Wf + Ug using previous inequalities we show that 

l|t//IL 




< 


'e-1 r 

4 Jo Jo 


dUg 


dz 


dzdt + 2 ll^s|li 2 (']o,T[;(L 2 (o,s))^ 


which leads to the announced property. ■ 

Theorem 16 For T small enough, 0) admits a solution in 

U/ (T\ - f ^ ^ I “ =v (z, t ), 

^ I V(;.,t)e]0,l[x]0,T[ 

Proof. Using (ITT|) with 4/ (ul) and T (Ug) implies 
||$(v[/K))-<F(4/(nf))||;.^(^) 

fT rl 

<T(aT + b) 



0 Jo 


Pf d 


r (1 — r^) dr 


dr 


dzdt 


Ko 
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because of m with a and b nonnegative. This leads to 




l|2 

W4T) 


<T{aT + b) ||T (nj) - T {u] 


(18) 


\\Wr{T) • 

Using (uni) with m we show that if T is small enough the mapping <l>oT is strictly contractant 
from hl4(T) to itself which prove the existence of a weak solution of ((3)) in the apropriate space. 


2.5 Back to the initial problem 


Let 


{L^ (0,1)) 


N 


du 


N 


dz 


gL2(0,1) , W,{T) = \ueLU0,T;W, 


fU (T) = I {u, v) e Wr (T) X W, (T) I u (1, z,t)=v (z, t), V (^, t) e [0,1] x [0, T]} , 
we have the 

Theorem 17 Assume 6is goes to 0 for i G {1,..., iV — 1}. The solution (Q/, Cis) of (jH) weakly 
converge in W (T) to the solution (^Cif,Cis^ of ([T]l. 

Proof. Take (p e W (T ), using (jH) we have 


T rl 



T rl 


Cif (r, 1, t) if (r, 1, f) r (1 — r^) drdt 



Cio (r) if (r, 0, f) r (1 — r^) drdt 



Cif^-r {1 — r'^) drdzdt + Pif f j f ^^^^-rdrdzdt 
lo Jo_^ Jo ^ Jo Jo Jo 

+ —[ Cis{z,T)(p{l,z,T)dz - — [ Ciso {z)ip{l,z,0)dz 

lis Jo Jis Jo 

lisJo Jo ot -fis Jo Jo dz dz 

A/ 


-/ ^i{Cis, ■ ■ ■ dzdt = 0. 

lis Jo Jo 

All the terms are bounded independently of 9is (cf. Lemma |2I). To pass to the limit in the 
non-linear term with r*, we use the fact that Cis is bounded in Lf ^ 0 , T; (L^ ( 0 , 1 ))^^ and 

the Theorem 2.1 of Krasnoselskii [8, page 22]. This allows us to let dis going to 0 for i G 
{1,...,A^ — 1} in the precedent variationnal formulation to obtain 


T rl 



T rl 


Cif (r, 1, t) ip (r, 1, f) r (1 — r^) drdt — 



Cio (r) ip (r, 0, t) r (1 — r^) drdt 



Cif-^r {1 — r"^) drdzdt + Pif f j f -—^^rdrdzdt 
'0 -'o ^ Jo Jo Jo dr dr 

+ —[ Cis{z,T)ip{l,z,T)dz - — [ Ciso {z)ip{l,z,0)dz 

Jis Jo Jis Jo _ 

-^r 0 d.it +4^ (i- 

lisJo Jo dt 'fis Jo Jo dz dz 

A/ 


-/ ^iiCis,---,CNs)T{^^^A)dzdt = 0. 

Jis Jo Jo 
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Making some particular choices for the test-function (p we deduce that \^Cif,Cisj is solution 
of 0 with the initial or boundary conditions ((21) • From (jS)) and the fact that the space of 
existence of Cif doesn’t depend on 9is we have that the solution (^Cif,Cis^ of Q is in PF (T) 
and the existence of a solution of (HD. ■ 


Remark 18 In the following we will write (Cj/, Cj^) instead of yCif,Cisj for the solution of 

(HD- 


3 Uniqueness of the solution 

Proposition 19 The system du admits an unique solution ^ m IF (T). 


Proof. Suppose there exists two solutions of problem ((T)), and set IF*/ = 

Cjj — Cfj, Wis = Clg — Cfg. For i G {1,iV}, (PF*/, PFs) is a weak solution of 


r dw,f 

dz 
dWi 


dt 


. ^ ^ 1 d ( dWif\ , 


{z,t) = 


z,t) 


(19) 


dr dz^ 

+6, (r, ...,Chs)- r. {Cl, C^)) (f t), 


with the initial or boundary conditions 

I Wif{r,0,t)= 0, PF/(l,2;,t)= Wis{z,t), ^^^{0,z,t)= 0, 

\0Ns^^{0,t)= 0, eNs^^{Ct)= 0, PF,,(^,0)= 0. 

Multiplying (fTTI|l i by r (1 — r^) PF/, integrating for r G ]0, T[, using and the fact that the 
r* are Lipschitz continuous with constant k we have 

{Wisf{z,T)dz<2kNj2 / {Wisfdzdt. 

i=i i=i 


But PFs(^, 0) = 0 so we use Gronwall’s lemma to conclude that PFs(2;,r) = 0 which implies 
that all the PFj/ and PFs are equal to 0 in their existence spaces. ■ 


4 Conclusion 

Starting from a non-stationary model of catalytic converter with cylindrical passage presenting 
many mathematical difficulties we have used some parabolic regularization technics and con¬ 
structed an appropriate functionnal space to prove that this problem admits one and only one 
solution for small time. 


12 













Acknowledgment. 

I would like to thank Professor Alain Brillard for giving me the subject of this paper for my 
PhD thesis, Professor Paul Deuring for his help in the proof of lemma 4, Professor Samir Akesbi 
for his advices and Professor Bernard Brighi for many useful discussions. 


References 

[1] BREZIS H., “Analyse fonctionnelle : theorie et applications.” Masson, Paris (1983). 

[2] CARASSO C., “Mathematical analysis of the model arising in study of chemical reactions 
in a catalytic cracking reactor.” Math. Comput. Modelling 18, No 2, 93-109 (1993). 

[3] CHIPOT M., “Elements of Nonlinear Analysis.” Birkhauser, Basel (2000). 

[4] DAUTRAY R. and LIONS J. L., “Analyse mathematique et calcul numerique.” Masson, 
Paris (1987). 

[5] HENRY D., “Geometric Theory of Semilinear Parabolic Equations.” Lect. Notes Math. 840. 
Springer Verlag, New York (1981). 

[6] HOERNEL J.-D., “Etudes theorique et numerique d’un modele non-stationnaire de catal- 
yseurs a passages cylindrique.” PhD thesis, Universite de Haute-Alsace (2002). 

[7] HOERNEL J.-D., “A non-stationary model for catalytic converters with cylindrical geom¬ 
etry,” Proceedings of the Fourth European Conference on Elliptic and Parabolic Problems - 
Rolduc and Gaeta 2001, World Scientihc, (2001), pp. 424-433. 

[8] KRASNOSELSKH M. A., “Topological methods in the theory of nonlinear integral equa¬ 
tions.” Pergamon Press, Oxford (1964). 

[9] KUFNER A., “Weighted Sobolev spaces.” Teubner-Texte zur Mathematik, 31, BSB B. G. 
Teubner Verlagsgesellschaft, Leipzig (1980). 

[10] OH S. H. and GAVENDISH J. G., “Transients of monolith catalytic converters. Response 
to step changes in feedstream temperatnre as related to controlling antomobile emissions.” 
Ind. Eng. Chem. Prod. Res. Devel. 21 29 (1982). 

[11] OH S. H. and GAVENDISH J. G., “Mathematical modeling of catalytic converter lightoff. 
Part HI ; Prediction of vehicle exhaust emissions and parametric analysis.” The American 
Institute of Chemical Engineers Journal, 31(6), (1985), pp. 943-949. 

[12] RYAN M. J., BEGKE E. R., and ZYGOURAKIS K., “Light-off performance of catalytic 
converters. The effect of heat/mass transfer characteristics.” SAE 910610 (1991). 

[13] ZEIDLER E., “Nonlinear fnnctional analysis and its applications.” Springer-Verlag, Berlin 
(1985). 


13 



